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Abstract 

"Top-down" models explain the observation of ultra high energy cosmic rays (UHECR; 
E > 5 ■ 10 19 eV) through the decay of very massive, long-lived U X particles" . If super- 
particles with masses near a TeV exist, X decays also lead to a significant flux of very 
energetic neutralinos, assumed to be the (stable or long-lived) lightest superparticles. 
There is a range of energies where neutrinos get absorbed in the Earth, but neutralinos 
can still traverse it. These neutralinos could in principle be detected. We calculate the 
detection rate in planned experiments such as OWL and EUSO. For bino-like neutralinos, 
which have been considered previously, we find detection rates below 1 event per Teraton 
of target and year in all cases; often the rates are much smaller. In contrast, if the neu- 
tralino is higgsino-like, more than ten events per year per Teraton might be observed, if 
the mass of the X particle is near its lower bound of ~ 10 12 GeV. 



1 Introduction 



The existence of ultra-high energy cosmic rays (UHECR), with E > 5 • 10 19 eV, remains a 
mystery pQ. They have been detected by every major cosmic ray experiment, but we do not 
know of any astronomical objects that has sufficiently strong electromagnetic fields extending 
over a sufficiently large volume to accelerate charged particles to the required energies. Nor do 
we understand how these particles, once created, can reach us, given their energy loss through 
scattering on the cosmic microwave background [2]. 

One radical idea [3] is that UHECR originate from the decay of very massive, yet long-lived 
X particles. Since one starts with very energetic particles, which lose energy first through 
parton showering and fragmentation, and later while propagating through the universe, these 
class of models are known as "top-down" models. The most energetic CR event that has been 
observed to date has E ~ 3 • 10 20 eV [3]. This implies a lower bound Mx ^ 10 12 GeV on the 
mass of the X particles. Since UHECR are observed today, the lifetime of X must be at least 
comparable to the age of the Universe. Several particle physics models containing candidates 
with sufficiently large mass and long lifetime have been suggested [HE]. Ways to produce these 
particles in the very early universe are discussed in [3J,[6]. 

Models of this type can be made compatible with all existing data, including the first data 
from the Pierre Auger observatory [7j. However, in order to decisively test these models, one 
has to find predictions that allow to discriminate between top-down and the more conventional 
bottom-up p] models. These two classes of models usually predict somewhat different spectra 
for photons and neutrinos at high energies, and/or different distributions of the arrival direc- 
tions. However, distinguishing between UHE photons and protons is nontrivial. Gigaton class 
neutrino telescopes now under construction should see some very energetic neutrinos if these 
models are correct [8] ; however, bottom-up models generically also lead to comparable neutrino 
fluxes. Anisotropies in the arrival direction can be expected [9], if X particles are distributed 
like (or even form the) Dark Matter in our galaxy; however, quantitative details depend on the 
distribution of matter near the galactic center, which is not well understood. 

These difficulties motivate the analysis of signals where bottom-up and top-down models 
make qualitatively different predictions. This may be possible if we postulate the existence of 
superparticles [10] at or near the electroweak energy scale. This assumption is quite natural in 
the given context, since supersymmetry is the only known way to stabilize the large hierarchy 
between Mx and the electroweak scale against radiative corrections^] Since Mx is much larger 
than the sparticle mass scale, X decays will produce large number of superparticles. This is 
true even if the primary decay of X only involves Standard Model (SM) particles; in this case 
superparticles will be produced in the subsequent parton shower [TH [12]. All these super- 
particles will decay into lightest superparticles (LSPs), assumed to be the lightest neutralino. 
In contrast, bottom-up models will produce a miniscule flux of superparticles. The reason is 
that the vast majority of UHE proton or photon interactions with matter only produces addi- 
tional light particles (in particular, light mesons and baryons); the cross section for producing 
superparticles remains very small even at these energies. 

This raises the question how one might observe these very energetic neutralinos. The crucial 
observation [T3] is that there is a range of energies where neutrinos get absorbed in the Earth, 

1 Note that "large" extra dimensions do not help here, since by construction the "fundamental scale" must be 
at least Mx in order to explain the observed UHECR; this is independent of the dimensionality of spacetime. 



1 



whereas neutralinos can traverse it with little or no loss of energy. The reason for this difference 
is the smaller neutralino-nucleon scattering cross section, and/or the smaller neutralino energy 
loss per interaction [2]. Note that neutralino interactions always lead to a superparticle in the 
final state, which will decay back into a neutralino. An interaction will therefore not change 
the total neutralino flux, but will shift it to lower energies, where it is (even) more difficult to 
detect. 

In this article we provide a detailed calculation of the neutralino event rates that one might 
expect in future cosmic ray detectors with very large target volumes, like OWL |15j and EUSO 
[16J. We improve on existing analyses [13J, [TTl [18] in several ways. We use neutralino spectra 
impinging on Earth calculated with the most complete code for X particle decays [12], where 
we analyze several different primary decay modes. We also carefully include the effects of 
neutralino propagation through the Earth, using the results of [H]. Our calculation of the 
event rates includes a cut on the visible energy deposited by a neutralino interaction; since this 
interaction again produces an invisible neutralino, the visible energy is usually significantly 
smaller than the energy of the incoming neutralino. Moreover, we investigate both bino- and 
higgsino-like neutralinos; the cross sections for the latter have also been computed in [H] . We 
find that higgsino-like neutralinos would in fact be much easier to detect; bino-like neutralinos 
most likely remain out of reach even for the planned EUSO and OWL missions. Finally, we 
calculate the neutrino background from the same model of X decays as the signal. 

The remainder of this article is organized as follows. The calculation of the event rates is 
described in Sec. 2. In Sec. 3 we present numerical results, and Sec. 4 is devoted to a brief 
summary and some conclusions. 



2 Calculation of Event Rates 

Neutralinos are produced along with protons, photons, electrons and neutrinos at the location 
of X decays, following a prolonged parton shower [TTl [T2] . We fix the normalization through 
the proton flux at 10 20 eV, which we take to be 

E 3 F p (E) = 1.6 • 10 24 eV 2 m- 2 s- 1 sr- 1 (1) 

at E — 10 20 eV. This roughly corresponds to the flux observed by the HiReS experiment 
[20J, which is somewhat smaller than that observed by AG AS A |19| . Note, however, that we 
ignore the contribution of photons to the UHECR flux. This is phenomenologically motivated 
by the observation that UHECR events seem to be proton-like, rather than photon-like [21] . 
Normalizing to the sum of the proton and photon fluxes would obviously reduce the predicted 
neutralino flux, and hence the event rate; depending on the X decay model, the reduction 
factor would roughly lie between two and five. On the other hand, we ignore all propagation 
effects. If most X decays occur at significant distance from our galaxy, which may well be true 
if X particles are confined to topological defects, both the proton and photon fluxes might be 
depleted by propagation, while leaving the neutralino (and neutrino) flux essentially unchanged. 
The presence of significant propagation effects would therefore increase the predicted neutralino 
flux on Earth. 

Neutralinos can interact with nucleons either through the exchange of a squark in the 
s— channel, or through the exchange of a Z° or gauge boson in the t— channel. In the 
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following we treat these two contributions, which essentially do not interfere |14j . in turn, 
before discussing the calculation of the neutrino-induced background. 

As explained in [T7J [TBI El) the s— channel contribution is dominated by the exchange of 
on-shell squarks. The event rate is given by: 

f^rnax /"X max fflm» q 1 T? dfjj^l y) 

^=E[ dx L 'y*?hp x) dy v - (2) 



Here, F^o is the differential neutralino flux, which depends on the neutralino energy as well 
as the matter deptfH X. The sum runs about all quark flavors q, and the first integration is 
over the visible energy E vis = E^a — E^o = yE^o . The factor 1/y appears because we 
integrate over the visible, rather than total, energy. The lower limit E m \ n on E W - 1S is determined 
by the energy sensitivity of the experiment, whereas the upper limit E m&K is determined by 
kinematics, E max ~ Mx/2; however, after propagation through the Earth the neutralino flux 
at the highest kinematically allowed energy is very small. The lower bound on the column 
depth, X min = 0.13 • 10 6 GeV 3 , corresponds to an angular cut of about 5% on the signal, 
i.e. we only count events that emerge at least five degrees below the horizon; this cut greatly 
reduces the neutrino background. X max = 2.398 • 10 6 GeV 3 is the maximal earth column depth, 
corresponding to neutralinos that emerge vertically out of the Earth. The kinematic maximum 
of the scaling variable y, for 2-body decays q — ■> q + Xi> is Z/max q = 1 ~ m \o/fn\- Since the 
maximal neutralino energy is finite, there should strictly speaking also be a non-vanishing lower 
bound on y\ note that we need the neutralino flux at E^o = E v ± s /y. An explicit expression for 
the differential cross section da s /dy can be found in [T3]. Finally, the constant factor V is given 
by 

V = 2nV cS e DC tN APw J D . (3) 

Here, V e s is the water equivalent (w.e.) effective volume, enc is the duty cycle (the fraction of 
time where the experiment can observe events), t is the observation time, Na = 6.022 x 10 23 g _1 
is Avogadro's number, p w = 10 6 gm~ 3 is the density of water, and Jd =\ dcos9/dX | is the 
Jacobian for the transformation cos8 — > X(cos6). 

The t— channel exchange diagrams predominantly lead to the production of heavier neu- 
tralinos or charginos in the final state [H], which we collectively denote by Xout- The visible 
energy therefore also depends on the Xout decay kinematics. The event rate can be written as: 

N t = / dE vis / dX / dy- J P-o(_^, X) (G*f(E Yis , y) + G?o(E via , y)j V . (4) 

^min J ^min 

y ~1 y V Al / 

Here we have written the contributions from charged and neutral currents separately. Each 
term is given by a convolution of a differential cross section for the production of Xout with 
the Xout decay spectrum. These convolutions are more easily written in terms of the variable 
z = E^o I E^q — 1 — y: 



*i L z\ az\ T dz 2 

(5) 



2 Matter depth X is costumarily given as a column depth, measured in g/cm 2 or, in natural units, in GeV 3 ; 
for the Earth, X e [0, 2.398 ■ 10 6 GeV 3 ] [22]. 
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Here z% = E^ ont /E^o describes the energy transfer from the incoming lightest neutralino to 
the heavier neutralino or chargino, and z 2 = E^o t /E^ out describes the energy transfer from 
this heavier neutralino or chargino to the lightest neutralino produced in its decay. z 2 is chosen 
such that z = Z\z 2 = 1 — y. Explicit expressions for the differential cross sections, and for the 
limits 2 m in,max, <2i,max in Eq.Q, can again be found in jT4j JU In principle one would need to 
include sums over Xont in Eq.([5]). In practice, however, a single neutralino and a single chargino 
dominate neutral and charged current reactions, respectively |14j . 

The event rates ([2]) and (@J depend on the neutralino flux after propagation through the 
Earth. Of course, the propagation effects also depend on whether s— or t— channel exchange is 
dominant. We treat these effects through straightforward numerical integration of the transport 
equations, as described in ffj] . 

The background is dominated by v T scattering through t— channel exchange of W or Z 
bosons. At the relevant energies electron and muon neutrinos get absorbed efficiently in the 
Earth. However, since v T interactions regenerate another u T , albeit at lower energy, r neutrinos 
can always traverse the Earth, although their energy may be reduced drastically. Again treating 
charged and neutral current processes separately, the background rate can be written as 

K = / dE vis / dX / dy~F v (^™, X) U y ,yj + NS°(E*, y) ) V , (6) 



where y — 1 — E Uiin / E UfiUt . In the case of NC scattering (Z— exchange) the entire visible energy 
results from the hadronic vertex. In case of CC scattering {W— exchange) we add the visible 
energy released in r decay to that produced at the hadronic vertex: 



N^ c (E vis ,y) = / — y - —6(z-z min )6(z max -z) 

' z\ dz\ T dz 2 



--l-y 

(7) 



This expression is formally very similar to Eq.(j5]), which also includes contributions to the 
visible energy from the decay of an unstable particle. This treatment is conservative since it 
ignores the fact that a r produced inside the target volume may decay outside of it. Moreover, 
if t production and decay both occur inside the target volume, it may be possible to use this 
"double bang" signature to remove these background events. On the other hand, we ignore the 
background from rs produced outside the target which decay inside the observed volume. This 
contribution should be smaller, since one would need higher neutrino energy to produce a given 
visible energy in this manner. Note that at the energies in question, r energy losses in rock or 
water are no longer negligible; this reduces the energy released in r decay even further. Recall 
that after propagation through the earth the v T flux is a steeply falling function of energy. 

The background rate ([6]) is proportional to the tau neutrino flux F v emerging from the 
Earth. The v T flux at the location of X decay is usually quite small [12J. However, due to 
near-maximal neutrino flavor mixing, the three neutrino fluxes impinging on Earth are very 
nearly equal, i.e. we take one third of the total neutrino flux, normalized according to Eq.(TjQ), 
as estimate of the incoming v T flux. 



,., 7VC CC NC CC i 

6 Note that the Gio ' of Eq.© are the integration kernels K~ ' of ref. [14], multiplied with the total 
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cross section for t— channel scattering. 
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As mentioned above, tau neutrinos may lose much of their energy while traversing the 
Earth. We solve the corresponding transport equations using the methods of ref. [14J. Since 
we are interested in very high energies, the tau leptons produced in CC v T reactions may 
lose a significant fraction of their energy before decaying. We therefore modified the standard 
treatment [22] in order to at least crudely estimate the effects of r energy loss in matter. We 
do this by formally treating this energy loss as additional scattering. To this end, we modify 
the integration kernel in the transport equation for v T as follows: 



da(E y ,z) f 1 da{E y ,zx)ldL{z x E v ,E"\ dEll 



a(E y ) dz J a(E y ) dz x L dE" 



(8) 

=E"/E 



Here E y — E/(l — y) is the energy of the incident neutrino that gives rise to a neutrino with 
energy E after the scattering, and the function dL(E Tjin , E TfiUt ) / dE T OUt describes the r energy 
loss. We make the very simple ansatz [23J 

dE 

—Z. = — (3pE T with (3 = 0.85 • lO^cmV 1 = const. (9) 

dz 

This implies E T (z) = E T (0)e~^ pz . We assume that all rs decay after traveling a distance 
Zdec = E T cr T /m T , where r r is the lifetime of the r lepton and c is the speed of light. Note that 
we estimate the average decay length from the r energy after propagation. This underestimates 
the decay length, and hence the effect of r energy loss. On the other hand, for E Ut < 10 10 GeV 
the ansatz (Q overestimates the energy loss [23J. Our approximation of a fixed decay length 
leads to 

dL ^' E P = 5 ( E " - E'eM-xE")) , (10) 

with constant k = (3pcT T /m T . The integral over dL/dE", which appears in Eq. ([8"]) , is then given 
by: 

where in the last expression E" has to be interpreted as a function of E', as determined by the 
argument of the 5— function. We can then evaluate the integral in Eq.([8~]): 



a(E y ) dz v ' -^y'--y-^y'^ Ey) dz > 



(12) 

z'=zi exp(KZiEy) 



The obvious advantage of our simplified treatment is that it does not necessitate the nu- 
merical evaluation of additional integrals. This would have been very costly, since the length 
scales involved in r energy loss and decay (a few km for E T ~ 10 8 GeV) are very much shorter 
than the v T interaction length in rock (~ 10 3 km for E Vt = 10 8 GeV) [23J. A more accurate 
treatment would therefore have required to use many more steps in X when integrating the 
transport equation; even with out simple treatment, or indeed without including the effects of 
t energy loss, calculating the v T flux emerging from Earth takes up to several CPU days. On 
the other hand, our simplified treatment can only give us an indication of the size of effects 
due to r energy losses. We find that the effect on the v T flux emerging from Earth is essen- 
tially negligible for E Vt < 10 7 GeV. This is also true for X > 0.3V max , since then the flux at 
E Ut > 10 7 GeV is negligible even if the r energy loss is ignored. However, it can reduce the v T 
flux by a factor of two or more at large E Vt and small X. 
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3 Results 



We are now ready to present numerical results. Earlier estimates [131 E] have shown that one 
will need at least teraton scale targets in order to detect hadronic interactions of neutralinos in 
top-down models. Currently the only technology that might allow to monitor such large targets 
is optical observation from space [T5l [16]. Here one detects the light, either from Cerenkov 
radiation or from fluorescence, emitted by very energetic showers in the atmosphere. The 
target is therefore quite thin: the neutralinos would have to interact either in the atmosphere 
itself, or just below it. One usually estimates an effective target thickness of 10 to 20 m w.e.. 
A teraton target then results if one can monitor O(10 6 ) km 2 simultaneously, which might be 
possible pj)i HE]. One drawback of this approach is that observations of this kind are only 
feasible on clear, moonless nights, leading to a duty cycle eoc in Eq-© of only about 10%. 
In our numerical results we therefore take a target mass of lTt, e^c — 0.1, and assume an 
observation time of one year. 

As shown in [13], the expected neutralino flux depends quite strongly on Mx as well as on 
the dominant X decay mode. Top-down models predict rather hard spectra, i.e. E 3 times the 
flux increases with energy. Fixing the (proton) flux at E — 10 20 eV therefore leads to smaller 
fluxes at E < 10 20 eV as Mx is increased. Moreover, if Mx is not far from its lower bound 
of ~ 10 12 GeV, much of the relevant neutralino flux is produced early in the parton cascade 
triggered by X decay, which is quite sensitive to the primary X decay mode. In contrast, if 
Mx 3> 10 12 GeV, in the relevant energy range most LSPs originate quite late in the cascade; in 
that case the LSP spectrum is largely determined by the dynamics of the cascade itself, which 
only depends on Standard Model interactions, and is not very sensitive to the primary X decay 
mode(s). 

Following ref. [13] we therefore study scenarios with Mx = 10 12 and 10 16 GeV, for four 
different primary X decay modes. In contrast to previous analyses [HI [TT], [18] we calculate 
the event rates for both bino-like and higgsino-like neutralinos. As explained in ref. [2] the 
former interact with hadronic matter almost exclusively through s— channel scattering, while 
the latter dominantly interact through t— channel diagrams. 

Finally, we present results for two different values of the minimal visible energy E m i n . Events 
with visible energy as "low" as 10 6 GeV might be observable via the Cerenkov light emitted 
by particles in the atmosphere with velocities exceeding the speed of light in air. On the other 
hand, the fluorescence signal (observed e.g. by the HiReS experiment [20]) can probably only 
be seen for energies > 10 9 GeV. In all cases we require the event to come from an angle at 
least five degrees below the horizon. This greatly reduces the neutrino-induced background, as 
explained earlier. 

We present results for higgsino- and bino-like neutralinos in Tables 1 and 2, respectively. We 
saw in ref. [13] that the cross section for neutralino-nucleon scattering depends only weakly on 
details of the sparticle spectrum if \\ is higgsino-like. In Table 1 we therefore only show results 
for one scenario with higgsino-like LSP. It has an LSP mass of 300 GeV, with the second 
neutralino and first chargino, which are produced predominantly in NC and CC scattering 
respectively, having masses of 310 and 303 GeV, respectively; the near-degeneracy of these 
three states is a consequence of these states all being higgsino-like, which in turn follows if the 
LSP is a rather pure higgsino state. 

As expected, we find much higher event rates for Mx = 10 12 GeV than for Mx = 10 16 
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Event rates for higgsino— like x' 
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Table 1: Predicted events rates per teraton and year (with duty cycle e^c — 0.1) for the 
scenario H2 of [H], where x? is higgsino-like, and for the v T induced background. Both signal 
and background depend on the mass Mx of the progenitor particle, as well as on the primary 
X decay mode. We show results for X decays into a first generation quark antiquark pair 
("gg"), into a first generation quark squark pair ("gg"), into a first generation lepton slepton 
pair ("//"), and into five quarks and five squarks ("5 x gg"). We only include events that emerge 
from an angle at least five degrees below the horizon. 



GeV. In the former case we also see that the predicted event rate depends significantly on the 
primary X decay mode, again as expected. The decay into a lepton plus a slepton turns out 
to be most favorable. The reason is that this decay mode leads to a rather small number of 
protons produced per X decay, or, put differently, to a large ratio of the LSP and proton fluxes 
[12J. Since we normalize to the proton flux, this then leads to a rather large LSP flux. This 
decay mode also leads to the hardest Xi spectrum. Since the primary X decay only involves 
weakly interacting (s)particles, parton showering carries away a relatively small fraction of the 
energy of the original particles. The original slepton will then eventually decay into a very 
energetic neutralino. As a result, increasing the cut on E vis by three orders of magnitude only 
reduces the predicted event rate by a factor of ~ 5 in this case. 

The second most favorable primary X decay mode is the one into five quarks and five 
squarks. Since we produce ten strongly interacting (s)particles already in the very first step, 
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each of which initiates an extended QCD shower, the final multiplicity is very large, but the 
fluxes are relatively soft. One then again needs a rather large normalization factor to reproduce 
the desired proton flux ([TJ at E = 10 11 GeV. Since the Xi spectrum is quite soft, increasing 
E min from 10 6 to 10 9 GeV now reduces the predicted signal by nearly two orders of magnitude. 

The worst case is X decay into SM quarks only. This gives a relatively hard proton spectrum. 
Moreover, superparticles are now only produced in the parton shower. This gives a small ratio 
of Xi to proton fluxes, and a relatively soft Xi spectrum. The fourth primary X decay we 
considered, into a quark and a squark, also leads to a relatively hard proton flux. However, 
since a superparticle is produced in the primary X decay, the Xi fhrx is larger, and significantly 
harder, than for X — > qq decays. 

We see that at least three of the four cases might lead to observable signals if Mx is near 
its lower bound, and if visible energies around 10 6 GeV can be detected. Of course, at that 
energy one expects a huge number of ordinary CR induced events, ~ 1 event per km 2 and 
second or (including the duty cycle) ~ 3 ■ 10 11 events per year in an experiment observing 
10 5 km 2 , as required for a teraton-scale target mass [24J. One will therefore need an excellent 
discrimination against such down-going events in order to extract the signal of at best a handful 
events per year. To that end one may need to sharpen the angular cut somewhat. This may 
also be desired to further reduce the v T induced background, which in this case is within an 
order of magnitude of the signal. Fig. 1 shows that for E m i n = 10 6 GeV, imposing a stronger 
angular cut will not reduce the signal very much. This is in accord with the results of ref . [14] . 
which show large neutralino propagation effects only for LSP energies well beyond 10 7 GeV in 
this case. Note, however, that typically E vis < 0.1E^o in for higgsino-like neutralino. 

On the other hand, only the most favorable scenario remains observable if E m i n has to be 
increased to 10 9 GeV. On the positive side, the u T induced background is now at least three 
orders of magnitude smaller than the signal, illustrating that the Earth can indeed be used as 
a filter. This is fortunate, since Fig. 1 shows that now the angular cut can be sharpened only 
at the cost of a significant reduction of the signal. However, in most cases one would need 
tens of Tt-yr to see a convincing signal even for Mx = 10 12 GeV; for Mx = 10 16 GeV and 
-Emm — 10 9 GeV, one would need Pt-yr of target mass times observation time! This would 
require monitoring virtually the entire surface of the Earth. The neutralino flux from decays of 
such very heavy X particle would remain invisible to teraton scale detectors even for a threshold 
energy of 10 6 GeV. Note that in this case the predicted event rate is almost independent of 
the primary X decay mode. The reason is that now the entire relevant energy range satisfies 
x = 2E/Mx <C 1, where the spectrum is determined almost uniquely by the dynamics of the 
parton shower |l2j. 

Table 2 shows event rates for bino-like neutralino. In this case the scattering cross section 
depends strongly on the squark mass [25| (TTJ [18]. We therefore show results for three different 
scenarios introduced in ref. [13], with first generation squark masses near 370, 580 and 1,000 
GeV, respectively. We see that the event rate remains below one event per year and teraton in 
all cases. This result seems much less promising than that of earlier studies [T3l [TT] . However, 
our rates are actually comparable to those of ref. [17] , once the differences in treatment are taken 
into account. To begin with, we assume that the X particles are distributed like Dark Matter, 
i.e. clump in our galaxy. Assuming a uniform distribution throughout the universe, as done in 
ref. [T7], increases the neutralino flux by about one order of magnitude [T3]. The reason is that 
such a uniform distribution suppresses the proton flux due to the GZK effect. One therefore 
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Figure 1: Angular dependence of the signal from higgsino-like neutralinos from primary X — > 11 
decays, and of the v T induced background, for two different values of the lower limit on the 
visible energy. 



has to increase the normalization in order to match the observed flux. A more or less uniform 
distribution of X particles could be achieved only if they are bound to cosmological defects, 
which nowadays are quite tightly constrained by analyses of cosmic microwave background 
anisotropies [26]. Moreover, we quote events per year, whereas ref.JTTJ finds about five events 
per lifetime of the experiment, taken to be three years. Finally, ref. [17] applies a cut (of 10 9 
GeV) on the total energy of the incident neutralino, whereas our cut is on the visible energy. 

We note that for -E m i n = 10 6 GeV, the ten body decay mode and X — > 11 decays now 
generally lead to similar event rates. The reason is that very energetic bino-like neutralinos 
lose energy considerably faster than higgsino-like neutralinos do: for rather light squarks the 
cross sections are comparable, but the energy loss per scattering is much larger for bino-like 
states, which produce a squark with nig 3> m^o, than for higgsino-like states, which produce a 
heavier neutralino or chargino very close in mass to the LSP. The 5xqq decay mode has a larger 
flux of softer neutralinos, which suffers less from propagation effects; for bino-like neutralinos 
this largely compensates the reduction of the rate due to the fact that the cross section is 
smaller at smaller LSP energy. However, if E v i S > 10 9 GeV is required, even the relatively 
softer LSPs produced from the ten body decay mode will typically scatter several times before 
reaching the detector. X — > 11 decays are then again more favorable, due to its initially much 
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Event rates for bino— like Xi 
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N D 3 


(7(7 

yy 


0.0191 


0.0192 


0.0118 


qq 


0.0471 


0.0528 


0.0388 


11 


0.3560 


0.5376 


0.5543 


5 x qq 


0.4567 


0.4779 


0.3051 


E vis > 10 9 GeV, M x = 10 12 GeV 


N D i 
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qq 
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0.00070 


0.00143 


qq 
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0.00314 


0.00701 
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0.14800 


5 X qq 
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0.01982 


0.03967 


E vis > 10 6 GeV, M x = 10 16 GeV 
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(7(7 
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00070 


00077 
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0.00079 


0.00117 


0.00062 


5 x qq 


0.00113 


0.00122 


0.00088 


E vis > 10 9 GeV, M x = 10 16 GeV 


N D1 


N D 2 


N D 3 


qq 


0.000006 


0.000058 


0.000140 


qq 


0.000005 


0.000047 


0.000107 


11 


0.000015 


0.000149 


0.000175 


5 x qq 


0.000006 


0.000067 


0.000161 



Table 2: Predicted event rates for bino-like LSP, for the same combinations of E m - m , Mx 
and primary X decay mode as in Table 1. We show results for the three different mSUGRA 
scenarios of [H], with first generation squark masses of about 370 GeV (Dl), 580 GeV (D2) 
and 1,000 GeV (D3). The background is essentially the same as in Table 1. 



larger flux of very energetic neutralinos. 

This brings us to a feature of our treatment which enhances the event rate compared to 
the numbers of ref . [T7] . In that analysis all neutralinos were discarded that interact even once 
before reaching the detector. This is not necessary, since this interaction will again yield a 
neutralino (from the decay of the produced squark), with typically about half the energy of the 
original LSP. Fig. 2 shows that this regeneration effect also leads to a much milder dependence 
of the final event rate on the cross section, and hence on the squark mass, than found in ref. p2] . 
Increasing the squark mass reduces the cross section, and hence the event rate for given flux. 
However, it also reduces the effect of neutralino propagation through the Earth, i.e. it increases 
the flux. These two effects obviously tend to cancel. As a result the event rate as function of rriq 
shows a rather broad maximum, the location of which depends on the cut on E V i S . A lower E V i S 
means that softer neutralinos can contribute. Since the cross section increases with neutralino 
energy, softer neutralinos can tolerate lighter squarks before suffering significant propagation 
losses. As a result, at smaller E m i n the maximum rate occurs for smaller squark mass. This 
effect is less pronounced for primary X —*■ 11 decays, since in this case the incident neutralino 
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spectrum is in any case rather hard, even if no cut on E vis is applied. 




Figure 2: Expected event rate due to bino-like neutralinos as function of the first generation 
squark mass, for two different primary X decay modes and two choices of the minimal visible 
energy See the text for further details. 



4 Summary and Conclusions 

In this paper we have calculated signal rates for the detection of very energetic neutralinos, as 
predicted by "top-down" models for the observed cosmic ray events at the highest energies. 
We use up-to-date calculations of the neutralino flux generated at the location of the decay of 
the superheavy particles, and of the effects due to propagation of the neutralinos through the 
Earth. We also for the first time treat the case of higgsino-like neutralino. 

We conservatively assume that the progenitor U X particles" are distributed like Dark Mat- 
ter, in which case most sources are "local" , i.e. effects of propagation through the interstellar 
or intergalactic medium are negligible. We then find detectable event rates in teraton scale 
experiments with duty cycle of ~ 10%, typical for experiments based on optical methods, only 
if the following conditions are satisfied: the lightest neutralino must be a higgsino, rather than 
a bino; Mx must be rather close to its lower bound of ~ 10 12 GeV; and the experiment must 
either be able to detect upgoing events with visible energy not much above 10 6 GeV, or most 
X particles undergo two-body decays involving at least one slepton and no strongly interacting 
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(s)particle. The good news is that in all cases we studied the signal is at least several times 
larger than the v T induced background, computed in the same X decay model. If Mx is near 
10 16 GeV and the LSP is higgsino-like, or Mx ~ 10 12 GeV and the LSP is bino-like, one will 
need (9(100) Tt-yr to collect a respectable event rate. In the worst case, with a bino-like LSP, 
Mx ~ 10 16 GeV and a threshold of the visible energy near 10 9 GeV, one would observe less 
than one event per year even if one monitored the entire surface of the Earth! These numbers 
improve by about one order of magnitude if X particles are distributed more or less uniformly 
throughout the universe; this might be expected if they are confined to cosmic strings or simi- 
lar topological defects. Recall, however, that scenarios with cosmic strings are constrained by 
observations of cosmic microwave anisotropies. 

These numbers only include interactions of neutralinos with nuclei. It has been claimed 
in Ref.|18] that bino-like LSPs should lead to a detectable signal in Gt class experiments 
(like IceCube [27]) through resonant production of sleptons. However, they estimate the rates 
assuming a neutralino flux close to the upper bound on the neutrino flux; the kind of model we 
investigate here yields fluxes that are several orders of magnitude smaller than this. Moreover, 
the visible energy in such events is relatively small, since only the decay of the produced slepton 
contributes. At the relevant energies the Earth does not filter tau neutrinos very well; so even 
if one concentrates on upgoing events, the background in potentially realistic X decay models 
is several orders of magnitude larger than the signal. 

Our overall conclusion is that next generation experiments, with effective target masses 
in the Tt range, would have to be lucky to observe a signal from neutralinos of "top-down" 
origin. Experiments with a relatively low energy threshold would stand a much better chance 
than those with high threshold. Unfortunately there are many reasonable X decay scenarios 
where the neutralino flux will remain invisible to such experiments. The goal of finding an 
experimentum crucis for top-down models may therefore remain elusive. 
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